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SUMMARY 

An  integral-transf orra  method  is  used  to  solve  the  elastic  problem  of  two 
col linear  cracks  at  the  edge  of  a  hole  in  a  sheet.  Both  uniform  loadings  on  the 
sheet  remote  from  the  crack  and  loadings  on  the  perimeter  of  the  hole  are 
considered.  Since  in  all  cases  the  loadings  are  symmetrical  about  the  crack¬ 
line,  only  the  opening-mode  stress  intensity  factors  are  non-zero;  these  factors 
are  calculated.  An  approximate  procedure  is  examined  for  obtaining  stress 
intensity  factors  for  two  cracks  at  a  hole  from  that  for  one  crack. 
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1  INTRODUCTION 

Cracks  in  the  engineering  materials  used  in  aircraft  structures  may  be 
present  at  manufacture  or  they  may  appear  during  the  service  life.  Indeed, 
recent  design  requirements  for  'damage-tolerant'  airframe  structures  are  often  | 

based  on  the  assumption  that  cracks  are  present  at  the  start  of  service  life. 

In  order  to  ensure  safety  and  to  optimize  inspection  schedules,  it  is  necessary 
to  know  both  the  residual  strength  of  the  cracked  structure  and  the  rate  at 
which  the  crack  will  grow  under  the  action  of  service  loads.  Both  of  these  depend 
upon  the  stress  intensity  factor  which  governs  the  stress-field  at  the  tip  of  the 
crack. 

Examination  of  failures  (see  for  instance  Kirkby*)  shows  that  both 
service  and  test  failures  frequently  occur  as  a  result  of  cracks  which  originate 
in  regions  of  high  stress  concentration  such  as  the  edges  of  holes  or  cut-outs. 

The  first  stress  intensity  factor  solution  available  for  radial  cracks  from  a 

2 

circular  hole  was  that  due  to  Bowie  .  His  solution  which  is  restricted  to  one 

crack  or  two  cracks  of  equal  length  is  not  very  accurate  for  short  cracks  owing 

to  the  limitations  of  the  method  used.  Knowledge  of  the  growth  of  short  cracks 

is,  however,  of  technological  importance,  since  much  of  the  fatigue  life  of  a 

structure  is  spent  while  the  cracks  are  short.  Recently  an  accurate  method  of 

solution  for  one  crack  at  the  edge  of  a  hole  has  been  obtained  by  Tweed  and 
3 

Rooke  ;  the  method  can  be  applied  to  any  configuration  with  a  stress  distribution 
which  is  symmetric  about  the  crack  line.  This  method  has  now  been  extended  by 

4 

Tweed  and  Rooke  to  problems  with  two  cracks  of  different  lengths,  a  common 
practical  configuration. 

In  this  Report  the  opening-mode  stress  intensity  factors  are  calculated  for 
two  collinear  radial  cracks  of  different  lengths  at  the  edge  of  a  circular  hole 
in  a  large  sheet  which  is  subjected  to  different  load  distributions.  Four  dis¬ 
tributions  are  considered;  they  are  (see  Fig  1) 

(i)  a  uniform  uniaxial  tensile  stress  remote  from  the  cracks  acting  in  a 
direction  perpendicular  to  the  crack  line, 

(ii)  a  uniform  biaxial  tensile  stress  remote  from  the  cracks  with  com¬ 
ponents  acting  parallel  and  perpendicular  to  the  crack-line, 

(iii)  two  opposing  forces  acting  perpendicular  to  the  crack  line  at  points 
on  the  edge  of  the  hole,  and 

(iv)  a  uniform  pressure  acting  around  the  perimeter  of  the  hole  (but  nr: 
on  the  crack  faces);  the  stress  intensity  factor  for  a  uniform 
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pressure  acting  around  the  perimeter  of  the  hole  and  on  the  crack 
faces  is  the  same  as  that  for  case  (ii)  if  the  two  components  of  the 
biaxial  stress  are  equal  to  the  pressure. 


In  section  6  it  is  shown  that  for  all  cases  the  predominant  parameter  in 
determining  the  stress  intensity  factor  is  the  total  crack-length,  ie  the  tip-to- 
tip  distance.  This  leads  to  a  relatively  simple  graphical  relationship  between 
results  for  two  cracks  and  results  for  one  crack.  An  approximation  which  has 
been  suggested  (see  for  example  Broek^)  for  obtaining  stress  intensity  factors 
for  two  equal-length  cracks  from  that  for  one  crack  is  examined  in  the  light 
of  this  relationship  and  its  accuracy  assessed. 

2  BASIC  THEORY 

The  hole  is  defined  in  plane  polar  co-ordinates  (p,  0)  by  0  p  £  R  , 

0  0  2ir  and  the  cracks  by  R  .<  p  $  cR  ,  0  =  0  and  R  £  p  $  bR  ,  0  =  ir 

(see  Fig  2).  The  crack  lengths  are  denoted  by  £j  and  £^  respectively  where 

X,  =  (b  -  1  )R  and  £  =  (c  -  1)R  .  In  general,  the  normal  stresses  along  the 

1  4 

crack  site,  in  the  uncracked  body  due  to  the  applied  loads,  are  given  by 


and 


a00(p’O)  =  '  P0f(r)  ’ 


aee<p’*)  =  "  po8(r)  * 


1  £  r  s  c 


I  <  r  ,<  b  , 


0) 


where  r  =  p/R  and  pn  has  the  dimensions  of  stress.  The  opening-mode 

,  .  u  .  4 

stress  intensity  factors  for  the  two  tips  at  (bR,  if)  and  (cR,  0)  are  given  , 
for  the  crack  of  length  £  ,  by 


K, 


n 


P0^ 


b  -  1 


h  (-b) 


and,  for  the  crack  of  length  £^  ,  by 


(2) 


/2 


p0^ 


c  -  1 


h(c)  . 


The  function  h(t)  must  satisfy  the  coupled  singular  integral  equations 
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-1 


h(t)M(r,t)dt 
/(-l  -  t)  (t  +  b)  17  /  /(c  -  t)  (t  -  1) 


+  i 


’  h ( t )M(r ,  t) 

•1  J  t  r  —  f-'Wt'  — 


t)  dt 


-b 

with  conditions 


f-g(-r) ,  -b  <  r  f  - 
l-f(r),  I  f  r  <  c  , 


h(  1 )  =  h(-l)  =  0  . 


The  kernel  M(r,t)  is  singular  and  given  by 

1 


M(r , t)  = 


t  -  r 


+  k(r,t)  . 


The  non-singular  part  k(r,t)  is  given  by 


(3) 


(4) 


(5) 


9  2  2  2 

k(r  ti  =  0  -  O  _  t(i  -  t  )  _  __t _  1 

3  2  1  -  rt  2 

t(l  -  rt)  (1  -  rt)z  r  t 


(6) 


3  NUMERICAL  PROCEDURE 

The  equations  (3)  and  (4)  are  reduced  to  a  system  of  simultaneous  linear 
algebraic  equations  by  using  a  procedure  developed  by  Erdogan  and  Gupta^.  Let 


|cos[(2j  -  l)ir/(2m)] 

j 

= 

1  j  2  j » » •  y  ui 

u.  =  * 

J 

u. 

J-m 

j 

m  +  1  , . .  • 

|Cos(jtr/m) 

j 

= 

1 , 2 ,  .  .  .  ,m 

V .  = 

J 

v.  , 

l  J"m 

j 

m  +  I  , . . . 

(b  -  l)u./2  -  (b  + 

0/2  , 

j 

= 

]  j  2  j  •  • « y  in 

t.  =  1 

J 

(c  -  1 )uj /2  +  (c  + 

0/2, 

j 

m  +  1 , . . . 

(b  -  l)v./2  -  (b  + 

0/2, 

j 

= 

1  1 2  y  •  *  • )  m 

r .  =  1 

'  J 

J 

(c  -  l)Vj/2  +  (c  + 

0/2, 

j 

= 

m  1  y  •  •  • 

(7) 


(8) 


(9) 


(10) 


Equations  (3)  and  (4)  can  now  be  reduced  to  the  following  linear  algebraic 
system: 


in 

o 


6 


2m 

Zh(t.)M(r.  ,t.) 

J  k  j 
j  =  l 


A  •  / 1  +  u* 

-  >  h(t.)(-l)J  - - J-  -  0 

®  Zj  j  V  1  -  u: 

J 


j=i 


2m 


-  £<rk)  , 


S  I  h(tj«<VtJ) 

j" 


:  1  h<ti n-niJr^  '  0  • 

j  =  l  J 


k  =  1 ,2, ...  ,m  -  1  ^ 


y  do 

k  =*  m  +  ],...,  2m  -  1 


J 


The  solution  h(t.)  of  these  equations  is  used  in  a  Gauss-Chebyshev  interpolation 
formula  to  determine  the  stress  intensity  factors  as  follows: 


for  the  crack  of  length  l ^  . 

4  LOADING  FUNCTIONS 

In  order  to  solve  equations  (II)  the  loading  functions  g(-x^)  and  f^x^) 
must  be  known.  These  are  derived  from  the  stress  fields  along  the  crack-site  in 
the  uncracked  configuration. 
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4. 1  Uniaxial  tensile  stress 

For  a  uniform  uniaxial  tensile  stress  a  acting  perpendicular  to  the 
crack-line  (see  Fig  la)  we  have7 


P0  =  a 


f(r)  =  g(-r)  =  1  +  — 

2r  2r 


4.2  Biaxial  tensile  stress 

For  a  uniform  biaxial  tensile  stress  of  a  and  ao  acting  respectively 
perpendicular  to  and  parallel  to  the  crack-line  (see  Fig  lb)  we  have7 


P0  =  o 


f(r)  »  g(-r)  =  f  j (r)  +  af2(r)  . 

The  function  ff(r)  is  given7  by  equation  (14)  and 


,  ,  ,  1  3 

^2^r  ~  2  ”  4 

2r  2r 


4.3  Point  forces 

For  two  opposing  forces  P  (force/unit  thickness)  acting,  diametrically 
opposite,  at  the  perimeter  of  the  hole  in  directions  perpendicular  to  the  crack- 

g 

line  (see  Fig  1c)  we  have 


>  =  -E- 

0  2R 


g(-r) 


"  r2  2  2 

L  o  +  J 


4.4  Constant  pressure 

For  a  constant  pressure  p  acting  on  the  perimeter  of  the  hole  (see 

g 

Fig  Id)  we  have 


8 


P0  =  P 


=  g(-r)  =  -=■ 


5  RESULTS 

Equations  (11)  have  been  solved  for  the  four  loading  systems  listed  in 
section  1  and  the  stress  intensity  factors  calculated  from  equations  (12). 

5. 1  Uniaxial  tensile  stress 

Opening-mode  stress  intensity  factors  have  been  obtained  for  the  crack  of 
length  l ^  (see  Fig  la)  as  a  function  of  for  various  values  of  £j/R  . 

The  results  are  tabulated  in  Table  1  and  shown  plotted  as  K^/ (a/nTI)  vs  X^/R  » 
for  various  values  of  £j/R  ,  in  Fig  2.  Table  1  also  contains  results  for  the 
special  case  when  the  two  cracks  are  equal  (£j  =  1^) .  The  limiting  value  of 
the  stress  intensity  factor  as  -+  0  is  given  by 


=  1.1215* 


where  *  is  the  stress  concentration  factor  at  the  edge  of  a  hole  diametrically 

opposite  a  radial  slit.  The  values  of  K  used  in  Table  1  were  taken  from 

.  9  11 

the  work  of  Wigglesworth  . 

5.2  Biaxial  tensile  stress 


The  opening-mode  stress  intensity  factor  for  the  crack  of  length 
(see  Fig  lb)  has  been  obtained  for  the  remote  biaxial  stress  field  of  a  per¬ 
pendicular  to  the  crack-line  and  a  parallel  to  the  crack-line.  The  results 
are  tabulated  in  Table  2  and  plotted  as  )  vs  X^/R  »  f°r  various  values 

of  £j/R  ,  in  Fig  3.  The  special  case  of  equal-length  cracks  (£^  = 
also  included  in  Table  2.  Stress  intensity  factors  K(a)  for  an  arbitrary 
biaxial  stress  field,  ie  a  perpendicular  to  the  crack  and  ao  parallel  to  the 
crack,  may  be  obtained  by  a  linear  combination  of  the  results  for  a  =  1  with  the 
uniaxial  results  (a  =  0)  as  follows: 


(1  -  a)K(O)  +  aK( 1 )  . 


(23) 
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5.3  Point  forces 

The  values  of  the  opening-mode  stress  intensity  factor  for  a  crack  of 
length  ^2  at  t^ie  edge  of  a  hole  subjected  to  opposing  point  forces  P  (see 
Fig  lc)  are  tabulated  in  Table  3  for  various  values  of  ;  results  for  the 
special  case  of  equal-length  cracks  (£j  =  Jt^)  are  also  included.  Curves  of 
K j/ (PQ  /n^)  VS  are  shown  in  Fig  4  for  various  values  of  £j/R  ;  in  this 

case  Pq  =  P/(2R),  the  ’bearing  pressure'. 

5.4  Constant  pressure 

For  a  case  of  a  constant  pressure  p^  in  the  hole,  but  not  in  the  cracks, 
the  stress  intensity  factor  has  been  evaluated  for  the  crack  of  length  i  . 
Values  are  tabulated  in  Table  4  and  plotted  as  K/  (p  /iri,  )  Vs  £.  /R  ,  in  Fig  5, 
for  various  values  of  JL^/R  .  Results  for  the  special  case  of  £  =  l  are  also 

included  in  Table  4. 

6  DISCUSSION 

The  accuracy  of  the  numerical  procedure  adopted  in  section  3  to  solve  the 
integral  equations  is  good;  the  results  are  accurate  to  <0.1%  for  m  =  20.  An 
advantage  of  the  method  used  in  this  Report  is  that  accurate  results  are 
obtained  for  the  important  region  of  short  cracks  without  the  need  for  any 
special  procedures  as  are  required  in  some  other  numerical  methods. 

Stress  intensity  factors  for  values  of  £j/R  not  contained  in  this  Report 
for  both  uniaxial  and  biaxial  stresses  can  best  be  obtained  by  interpolation  of 
the  function  K/ (o/rra)  ,  where  2a  is  the  total  crack-length  from  tip-to-tip, 

ie 

2a  =  £  +  2R  +  l  .  (24) 


o 


Values  of  K^/(a/ira)  are  tabulated  as  a  function  of  l^/R  for  the  various 
values  of  Jl|/R  for  both  uniaxial  and  biaxial  applied  stress  in  Tables  5  and  6 
respectively.  It  can  be  seen  from  the  tables  that  for  a  given  value  of  i^/R 
the  value  of  Kj/(a/rra)  is  not  very  dependent  on  £^/R  .  This  means  that  the 
total  length  of  the  crack  2a  (which  includes  the  hole  diameter)  is  an  important 
parameter  in  determining  the  stress  intensity  factor.  The  stress  intensity 
factor  can  be  written  as 


K 


I 


Y 


n/TTs  +  v  )  n 


2 


1 


2 


(25) 
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for  two  cracks  and  as 


Kj.  =  Y1a/7r"(2R  +  1^/2 


for  one  crack.  Thus  from  equations  (25)  and  (26), 


Kj(two) 

Kj(one) 


la  IJEHU 

Yj  J  S. j  +  2R  +  i2  * 


The  approximate  relationship, 


K  (two)  - 


2R  +  £2 
'  +  2R  +  l2 


(one) 


(26) 


(27) 


(28) 


has  been  suggested  by  several  authors  for  the  case  £  =  £  (see  for  example 
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Broek  );  this  relationship  assumes  that  Y^Yj  *s  uu^Y*  From  the  values  of 

Y^  and  Yf  ,  given  in  Tables  5  and  6,  the  ratio  Y2/Yj  has  been  calculated  and 

is  plotted  in  Fig  6  for  both  uniaxial  and  biaxial  applied  stresses.  The  curves 

of  y2^Yi  ys  ^2^  are  Plotte^  ^or  various  values  of  £j/R  .  It  can  be  seen 

that  ^2^1  c^ose  to  unity  for  the  data  considered  and  therefore  equation  (28) 

can  be  applied  to  the  crack  of  length  £2  for  £j  =£  £2  *  For  t'ie  uniaxiai  stress 

case  the  deviation  from  unity  is  less  than  5%  for  £2/R  <  0.5  and  less  than 

10%  for  £,j/R  <  2.0.  For  biaxial  stressing  the  ratio  Y2^Y1  ^as  a  minimuin  at 

£2/R  “  0.5  ;  for  all  the  data  considered  the  deviation  from  unity  is  between 

+3%  and  -4% , 

The  Y-functions  defined  in  equations  (25)  and  (26)  relate  the  stress 
intensity  factors  for  a  crack  at  an  unloaded  hole  to  that  for  an  isolated  crack 
(ie  a/ira)  .  In  a  similar  manner  Y-functions  can  be  defined  for  a  crack  at  a 
loaded  hole.  Thus 


K 


I 


/  2 ,  +  R 

J  7^-r  (2a  =  *,  +  2R  +  V 


(29) 


where  P  is  the  force  per  unit  thickness, 
of  the  tips  of  an  isolated  crack  of  length 
on  opposite  faces  of  the  crack,  a  distance 
£j  +  R  from  the  other  tip,  is  given  by 


The  stress  intensity  factor  for  one 
2a  subjected  to  forces  P  acting 
£2  +  R  from  the  considered  tip  and 


1 05 


p 


1 1 


(30) 


The  values  of  the  Y-f unctions  are  tabulated,  as  a  function  of  for  various 

values  of  £  /R  ,  in  Tables  7  and  8  for  the  point  force  P  and  the  internal 
pressure  p^  respectively.  In  the  case  of  the  internal  pressure  the  force  is 
equal  to  2p^R  .  As  in  the  case  of  remote  stressing  conditions  the  Y-functions 
are  not  very  dependent  on  £j/R  ;  thus  these  data  can  be  used  to  interpolate 
values  of  the  stress  intensity  factor  for  any  value  of  £j/R  not  contained  in 
the  results  reported. 

From  equation  (29)  it  follows  that,  for  two  cracks. 


2(£  +  R) 


U,  +  2R  +  *  )(£  +  R) 


(31) 


and  for  one  crack 


2R 


(2R  +  *  )(fc2  +  R) 


(32) 


Thus 


K  (two)  Y  /O’  +  R)  (2R  +  £  ) 

I  _  Z  i  1 _ z_ 

K^XoneT  ~  Yj  V  R(«(  +  2R  +  J^) 


(33) 


Equation  (33)  is  similar  in  form  to  equation  (27)  which  suggests  a  similar  approxi¬ 
mation  to  that  in  equation  (28)  for  estimating  the  stress  intensity  factors  for 
the  crack  of  length  £ when  l ^  ^  0  from  that  for  a  crack  of  length  when 

=  0  .  The  approximation  is 


K  ( two)  /(«!  +  R)  ( 2R  +  t>.y) 

(orTeJ  ~  V  R  (9~  +  2R  + 


(34) 


This  approximation  assumes  that  Y.,/Y  —  1  ;  the  errors  involved  in  this 

assumption  can  be  seen  from  Fig  6  where  Y  /Y  vc  R  is  plotted  for  various 


values  of  £(/R  for  both  point  loads  and  uniform  pressure  on  the  hole.  For 
point  loads  the  maximum  deviation  of  Y2/Yj  ^rom  unity  occurs  for  Jlj/R  =  0.5; 
the  deviation  varies  from  -4%  for  small  i^/R  values  to  -7%  for  large  values 
of  •  For  uniform  pressure  the  deviations  are  between  ±6%  for  <  1.0 

and  £j/R  <10;  for  large  values  of  Jl^/R  deviations  are  greater,  -13%  for 
£  /R  =£2/R=  10  . 


The  simple  approximate  relationships  for  obtaining  opening-mode  stress 

intensity  factors  for  two  cracks  from  those  of  one  crack  for  both  remote  and 

local  loading  would  still  apply  if  the  loading  was  not  symmetrical  about  the 

crack-line.  Crack  problems  with  asymmetrical  loading  about  the  crack-line  are 

more  difficult  to  solve,  particularly  if  there  are  two  cracks  of  unequal  length. 

However,  some  solutions  for  a  single  crack  under  remote  loading*^’1*  and 

.  .  1112. 

localized  loading  *  are  available  and  could  be  used  to  obtain  approximate 
opening-mode  stress  intensity  factors  for  two-crack  problems. 


7  CONCLUSIONS 

(1)  Accurate  opening-mode  stress  intensity  factors  have  been  obtained  for 
cracks  at  the  edges  of  a  circular  hole  in  sheets  subjected  to  either  remote 
tensile  stresses  or  loads  on  the  perimeter  of  the  hole. 


(2)  An  approximate  technique  for  obtaining  stress  intensity  factors  for  two 
cracks  from  values  for  a  single  crack  has  been  investigated  and  the  errors  in 
the  approximation  shown  to  be  small  (a  few  per  cent). 

(3)  The  approximate  technique  could  be  extended  to  other  configurations  for 
which  only  single  crack  results  are  known. 
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crack-length  (2a  =  tip-to-tip  distance) 
crack  tips  located  at  r  =  b,  c 
loading  functions  (dimensionless) 
loading  function  (dimensionless) 
function  defined  by  equation  (3) 
end-point  values  of  h(t) 

suffices;  j  =  1,  2,  . m,  k  =  1,  2,  . m  -  1 
opening-mode  stress  intensity  factor 
values  of  at  tips  b  and  c  respectively 

stress  aoneentration  factor 

non-singular  function  defined  by  equation  (6) 
crack-length  =  R(b  -  1) 
crack-length  =  R(c  -  I) 

number  of  integration  points  (section  3) 
singular  kernel  in  equation  (3) 
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dimensionless  radial  coordinate  (=  p/R) 
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sidered.  Since  in  all  cases  the  loadings  are  symmetrical  about  the  crack-line, 
only  the  opening-mode  stress  intensity  factors  are  non-zero;  these  factors  are 
calculated.  An  approximate  procedure  is  examined  for  obtaining  stress  intensity 
factors  for  two  cracks  at  a  hole  from  that  for  one  crack. 


